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Abstract 

We consider the condensation energy, E c of strongly coupled, magnetically-mediated supercon- 
ductors within the context of the spin-fermion model. We argue that although experimentally 
obtained values of E c are of the same order of magnitude as would be expected from BCS theory 
in optimally and overdoped cuprates, this agreement is coincidental. The actual physics behind 
the condensation energy is much richer. In particular, we argue that it is vital to take both the 
fermionic and bosonic contributions to the condensation energy into account when considering such 
materials, and that it is only the sum of the two contributions, E c , which has physical meaning. 
Both the experimental and our theoretically calculated condensation energies exhibit a decrease 
with further underdoping past optimal doping, e.g. in the mid to strong coupling regime. Below 
optimal doping, the physics is qualitatively different from BCS theory as the gain in the conden- 
sation energy is a result of the feedback on spin excitations, while the fermionic contribution to 
E c is positive due to an "undressing" feedback on the fermions. We argue that the same feedback 
effect accounts for a gain in the kinetic energy at strong coupling. 

PACS numbers: 74.25.-q, 74.72.-h, 61.12.-q 
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I. INTRODUCTION 



Undertanding the origin of the condensation energy is an important step towards iden- 
tifying the mechanism of high temperature superconductivity in the cuprates. In a BCS 
superconductor, the condensation energy E c - the energy gain in a superconductor com- 
pared to the normal state at the same T, smoothly increases below T c and at T = 
reaches E^ cs = —VNfA 2 /2, where V is the volume, A is the superconducting gap and 
Nf = mpp I (2ir 2 h 3 ) is the fermionic density of states. [1] The decrease in the total energy 
upon pairing results from a fine competition between an increased kinetic energy and a de- 
creased potential energy, both of which are much larger than E c . The BCS condensation 
energy can be experimentally extracted from the jump of the specific heat at T c as within 
the BCS theory C s - C n « 6.08E C /T C . 

Since the fermionic density of states is only weakly dependent on doping, the application 
of the BCS formula for the condensation energy to cuprate superconductors would imply 
that E c and A 2 scale in the same way. However, the measured gap increases monotonically 
with decreasing doping [2], while the jump of the specific heat has a non-monotonic doping 
dependence. In the overdoped regime it initally increases with reduced doping, but below 
optimal doping further underdoping leads to a decrease in the specific heat jump [3]. This 
discrepancy between the trends in C s — C n and in A as functions of doping clearly makes 
the applicability of the BCS formula to the cuprates questionable. 

In view of this fact, and the general consensus that the cuprates are strongly coupled 
systems, it comes as no surprise that effort has already been made to explain the cuprate 
condensation energy as a result of non-BCS physics. Scalapino and White [4] conjectured 
that at strong coupling, the dominant contribution to the condensation energy comes from 
a feedback effect on the magnetic excitations of the system. On the other hand, Norman et 
al argued [5] that the condensation energy likely has an electronic origin, and is driven 
by a gain in the kinetic energy which at strong coupling is negative (in contrast to BCS 
theory) because of a strong "undressing" of fermions which bear a greater resemblance to 
free particles in the superconducting state than they do in the normal state. Similar ideas 
were expressed by Hirsh and Marsiglio [6]. A different idea, related to the lowering of the 
Coulomb energy in the superconducting state has been proposed by Leggett [7]. 

In this communication we argue that these apparently disparate viewpoints are in fact 
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consistent with each other, and describe the same strong coupling physics. We argue that at 
strong coupling, the relation between E c and A is qualitatively different from BCS theory and 
is consistent with the experimental trends in the underdoped cuprates. We show furthermore 
that the strong coupling effects are in large part the result of mutual feedback between the 
fermions and bosons. We make the case that the contributions to the condensation energy 
from these two channels may not be considered independently from each other, and that 
only the sum of the two contributions has physical meaning. 

Our point of departure is the general equation for the free energy of an interacting electron 
system, derived by Luttinger and Ward [8] for the normal state and extended to a super- 
conductor by Eliashberg [9]. We first review strongly coupled phonon superconductors. 
We discuss the validity of the Eliashberg formula and describe the BCS, Bardeen-Stephen 
[10] and Wada [11] expressions for the condensation energy (for a review, see Ref. [12]). 
The BCS expression neglects both bosonic and fermionic self-energies (apart from a trivial 
renormalization of the dispersion). The Bardeen-Steven and Wada expressions include the 
fermionic self-energy and its change between the normal and the superconducting states, but 
neglect the feedback from superconductivity on phonons. The two expressions are generally 
considered to be identical, however we argue that this is only true as long as the feedback 
from fermions onto phonons is neglected. We review Bardeen and Steven's arguments for 
the validity of this approximation for phonon superconductors. 

We next modify the Eliashberg equation for E c to the case of magnetically mediated 
superconductors. We argue that in qualitative distinction to the phonon case, the feedback 
effect from the superconductivity onto the bosons may not be neglected if the pairing is 
magnetic. We derive the relevant equations and explicitly compute E c assuming that the 
pairing is due to spin-fluctuation exchange and is described by the spin-fermion model. We 
show that when the feedback on bosons is non-negligible, the calculations require care as 
both electronic and spin parts of the condensation energy are ultraviolet divergent. We 
explicitly show that these divergencies are cancelled out between the two terms, and the 
total E c (which turns out to be the only physically meaningful quantity) is free from diver- 
gencies. We furthermore demonstrate that one can avoid the divergencies by performing the 
computations in real frequencies. We apply the results to the cuprates and show that our 
theoretical E c agrees with the data both in magnitude and in the doping dependence. We 
view this agreement as a support for the spin-fluctuation scenario for the cuprates. 
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II. PHONON SUPERCONDUCTORS 



The condensation energy E c is the difference between the free energies in the normal and 
superconducting states. 



In the Greens function technique one evaluates the grand free energy Vt. The difference 
F n — F s coincides with Vt s — Q n provided that the chemical potential \i does not change 
between the two states. In the Eliashberg theory that we will be using, this is the case as 
the Fermi energy is assumed to be much larger than A, and the corrections to /i due to 
pairing, which scale as powers of A/E F , are neglected [10]. 

There are two ways to compute E c using Green's functions. The first approach is to use 
the general formula for the ground state energy of the interacting fermionic system in terms 
of the integral over the running coupling constant [13]: 



where f2 is the Free energy for free fermions, £ is the full self-energy, and G = (i(u + 
E(fc, A')) — efc) -1 is the full Green's function. Both are functions of a 4 vector of momentum 
and Matsubara frequency k = (k, ou n ) and the running coupling constant A'. Note that here 
and below we define £ with an extra factor of %. 

Eq. (2) is applicable to both normal and superconducting states (for the latter, E(fc, A') 
has a pole), however, is not convenient for our purposes as numerical calculations then re- 
quire solving the Eliashberg equations for a large set of coupling values. Still, we tried this 
formalism, and will discuss the results in the Appendix A. The approach which we will be 
using in the bulk of the paper was suggested by Luttinger and Ward [8] who demonstrated 
that it is possible to re-express Q in the normal state via a series of closed linked skeleton di- 
agrams with fully dressed fermionic and bosonic propagators. Their approach was extended 
to the superconducting state by Eliashberg [9]. We refer the reader to Refs.[8] and [9] for 
the details of the derivation, and here just present the result. In the superconducting state, 
f2 has the form 



E = F - F 



(1) 





n = 



IT £{i log[4 + t\k) + <$>\k)\ - iT,(k)G(k) + 



k 



+ 



^{log^^ + n^z^)} 
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+ T 2 J2 al-k>{G(k)D(k - k')G(k') + F(k)D(k - k')F(k')} + ... (3) 

k,k' 

The last term in the above equation is the sum of the first two closed linked skeleton 
diagrams, the dots stand for higher order diagrams. The functions G(k) and F{k) are the 
real and anomalous Greens functions given by 

e k + iE(fc) 



G(k) 



4 + E(A;) 2 + <l> 2 (A;)' 



F(k) = i 39 (4) 

where <§>(k) is the pairing vertex and S(/c) = u n + E(/c). The conventionally defined pairing 
gap A{k) is the ratio of the anomalous vertex and the self-energy: A(k) = &(k)u m /Tl(k). 
Finally, a k is the electro n-phonon coupling, and D{q) is the dressed phonon propagator given 
by D^ 1 (q) = D 1 (q) — 11(g) where 11(g) is the polarization bubble, and D (q) is the bare 
propagator. A careful reader may note that we define D(q) as dimensionless quantity (see 
(3), while the actual phonon propagator has a dimension of inverse energy. This does not 
cause problems, however, as D appears only in a combination with a 2 , and the extra overall 
energy factor can be absorbed from a 2 .. We will be using Eq. 3 as the point of departure 
for our analysis. 

Eq. 3 is quite general. Its use for the calculation of E c takes into account not just the 
introduction of the pairing vertex $ but also changes in the fermionic self energy E and the 
polarization bubble II. Thus it calculates the contributions to E c from both the fermions 
and bosons. 

The approximations made in the phonon case are associated with the smallness of the 
sound velocity v s compared to the Fermi velocity vf- It turns out that higher order closed 
linked diagrams form series in powers of Xv s /vf, where A is the dimensional coupling constant 
which scales with a k - k > (see below). The approximation, attributed to Migdal [14] and 
Eliashberg [15], is to neglect 0(Xv s /v F ) terms without assuming that A is by itself small. 
This would imply neglecting all terms labled as dots in (3). 

The physics behind the Migdal-Eliashberg approximation is best revealed by analyzing 
the perturbation series for E(/c) and . These series can be separated into terms that scale 
as powers of A and terms that scale as powers of Xv s /vf- The perturbative series in powers of 
\v s /vp arises from the actual electron-phonon scattering, and the factor v s /vp results from 
the fact that in this process electrons are forced to vibrate at phonon frequencies far from 



their own resonance. This gives rise to vertex corrections (at typical pairing frequencies), 
and to an equal renormalization of fermionic u> and e^, i.e. the quasiparticle residue is 
renormalized, but the quasiparticle mass remains a bare one. The vertex correction diagrams 
for and $(&), embedded into the closed linked diagrams, gives rise to higher-order 

skeleton diagrams for the Free energy, and it is these diagrams which are dropped. 

The terms that form series in A are different and they are not considered to be small. 
In the normal state, they can be understood as coming from phonon-induced interactions 
between electrons and their own zero-sound collective modes. These terms do not contribute 
to the vertex renormalization (at typical frequencies for the pairing), and gives rise to a 
that only depends on frequency. The 0(A) terms also contribute to the pairing problem and 
in the superconducting state give rise to $(&) that again depends only on frequency. 

Neglecting the higher order skeleton diagrams, Eliashberg obtained the closed-form ex- 
pression for Q s given in Eq. 3. The closed form, coupled equations for the fermionic self- 
energy E(/c) = S(c<j) and the pairing vertex <£>(&) = <&(ui) then follow from the condition that 
il s given by (3) is stationary with respect to variations in E(fc) and $(&;). The conditions 
5Q s /5E(k) = 5n s /5<S>(k) = yield 

t(k) = to + iTj2at_ k ,G(k')D(k - k') 
k' 

<!>(k) = -iTj2a 2 k _ kl F(k')D(k-k') (5) 

k' 

The Migdal- Eliashberg approximation has two further implications. First, the fermionic 
dispersion may be approximated by e k = v F {k — k F ), as typical pairing frequencies should 
be much smaller than Ep if \vsjvp is to be small. Second, the momentum integration over 
k' can be factorized: the integration over momenta transverse to the Fermi surface involves 
only normal and anomalous fermionic propagators, via e k , while the integration along the 
Fermi surface involves only the bosonic propagator in which one can set \k\ — \k f \ — kp, i.e., 
deviations from the Fermi surface are neglected leading to a bosonic propagator dependant 
only on frequency. Corrections to this approximation again scale as Xv s /vp. Finally, it is 
assumed (without justification) that the Fermi surface is isotropic in the sense that vp is 
considered to be independent of k. Under these approximations, the momentum integration 
in the equations for and $(cj) can be performed exactly. Approximating the momentum 



6 



sum by an integral over energy and an associated density of states Nf as 

/oo 
de (6) 

, inserting equations 4 into equation 5 and performing the integration we obtain 

XV) 



= u + TrTN f '£al_ u ,D(uj-u') 



$M = nTNf^al^Diu-u;') ttL (7) 

'A 2 (a') + E 2 (u') 



An essential feature of the above equations is that apart from the assumption of the 
smallness of Xv s /vp, it is assumed assumed that the phonon polarization bubble H(k), which 
accounts for the effects of the electrons on phonons, may be neglected. Analagously to the 
derivation of equation 5, an expression for Tl(k) can be formally obtained from the feature 
that Q s given by (3) is also stationary with respect to variations in U(k). The condition 
5Q s /5U(k) = yields 

U(k) = -2TD(0)J2^l-AG(k)G(k - k!) + F(k)F(k - k')] (8) 

k' 

The fact that H(k) is irrelevant to the phonon problem is not immediately apparent and 
this issue must be considered carefully. In the normal state, 

n(Jfe) oc A— — ^f- (9) 
v F cu D |k| 

reflecting the decay of a low-energy bosonic mode into a fermionic particle-hole pair. This de- 
cay term obviously cannot be neglected at the lowest frequencies as it accounts for the leading 
low-frequency dependence of D{uj m ). I.e. for an Einstein phonon, D (u m ) = 2uj 2 d / {u^ + u^) 
where ujj is the Debye frequency (we recall that D{oj m ) is dimensionless, hence Z^~ 1 (cj m ) = 
D 1 (uj m ) — n(/c) ps 1 — n(fc) at the lowest frequencies. However, for A = 0(1), the frequen- 
cies relevant to the pairing problem are of order ujd, typical momenta are of order pp, and 
hence H(k) at typical frequencies is small to the same extent as Xv s /vf and other similar 
terms in the Eliashberg theory. In the superconducting state, the low-energy phonons are 
gapped, and at frequencies smaller than A, 

n(Jfe) oc A— ^ — jf-. (10) 
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This difference between H(k) in the normal and superconducting states reflects a fundamen- 
tal change in the bosonic dynamics at small frequencies due to the gapping of low-energy 
fermionic excitations. Still, however, at typical frequencies for the pairing, uj ~ A ~ u D , 
H(k) ~ Xv s /vp <C 1, and the polarization operator can be neglected compared to the bare 
phonon propagator. This was shown explicitly by Bardeen and Stephen [10] who evaluated 
the contribution to E c from the change of the bosonic dynamics between normal and super- 
conducting states and demonstrated that in the phonon problem the bosonic piece in E c is 
small compared to the electronic piece and hence can be safely neglected. 

We will see below that neglecting the bosonic contribution to E c is not a justifiable 
approximation if one is dealing with magnetic superconductors where the bosonic mode is a 
collective mode of the fermions. This is the main difference between phonon and magnetic 
superconductors, and we will examine this issue in the next section. 

For completeness, we present several useful forms for £l s and E c for the phonon case. For 
simplicity, we assume that the electron-phonon coupling is independent of frequency, and 
the phonon spectrum consists of a single Einstein boson with a frequency ujb ~ v s pf, i.e., 



ex./,) — cy, 



o .2 

D(q) = D{u m ) = -5-^- (11) 



Integrating over momentum in Eq. (3) as before, we obtain 
n = -N f {2irTj2 



+ tw v + 1 — } (12) 

— "2 — 1— $2 /g2 7 $2 [U m -U m >y +UJ Z D 



where we introduced a 2 = 2a 2 NfU> 2 :) . The dimensionless coupling A introduced above is 
related to a 2 as 



a 2 

X = ^ (13) 



Eq. (12) may be simplified by making the standard substitutions S Wm = u m + S Wm = 
u m Zu m and $ Wm = A Wm £ Wm /w m = A^ m Z^ m . Substituting these forms into (12) we obtain 



n = -N f {2nT^ 



U 2 m 



^ra 2 + Al 



+ T\ 2 a 2 £ - ^m' + A^A^, 1 



Thus we see that the free energy for a phonon superconductor is dependent only upon the 
form of the gap A Wm and is not explicitly upon the self energy E Wm . 

A comment is in order here. The Luttinger-Ward result for the Free energy as a series of 
closed linked skeleton diagrams is, strictly speaking, only valid for the minimum of Q, i.e., 
for the self-energy that satisfies the stationary condition. Otherwise, the Luttinger-Ward 
generating functional (which is what they actually calculated) does not necessarily coincide 
with the Free energy. In the normal state, this does not cause a problem with Eq. (14) as 
£l N does not explicitly depend on Sjy. (We use the capital subscripts "S" and "N" to denote 
"superconducting state" and "normal state" respectively. Lower case "n" (and "m") refer to 
Matsubara frequencies.) In the superconducting state, Eq (14) does not imply that Qs is a t 
a minimum, i.e., 5Qs/$A — 0. We didn't analyze in detail the corrections to Eq. (14) which 
would stem from the difference between the Luttinger-Ward-Eliashberg functional and the 
actual Free energy, but the estimates show that these corrections would be again be small 
in Xv s /v F . If this is the case, then Eq. 14 can be used for the study of the profile of the free 
energy, i.e. how it evolves for different solutions of the gap. In particular, we verified that 
for the BCS case, the expansion of Eq. 14 near A = to order A 2 yields the sign change of 
the slope at exactly the BCS transition temperature T c . 

For the remainder of this communication, we will only be considering Q evaluated at 
the equilibrium solution of A Wm when the applicability of the Luttinger-Ward formalism is 
rigorously justified. To this end, the Eliashberg equation for the equilibrium solution A^ 
may be obtained by minimizing equation 14 with respect to A. This yields 



This solution for A must then be substituted into equation 14 allowing the calculation of 
the free energy at its minimum, and hence E c = Q$ where fl^ is obtained from Q$ by 
setting A = 0. Performing the computations, we obtained Eq. (17) below. 

There is, however, a simpler way to proceed towards E c . In this approach, one specifies at 
the outset that one is only interested in the Free energy at equlibrium. In this case, the fact 
that equation 3 should be stationary with respect to variations in E and $ can be invoked 
even before the momentum integration is performed [10-12]. Substituting equation (5) into 



A, 




1 



1 



(15) 
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equation 3 and droping the phonon piece, we obtain 

A = -T^{log^ + S^ + ^J-^ m G(A;)+^ m F(A;)} (16) 
v 

Introducing a density of states Nf, integrating over and subtracting the normal state 
result from the superconducting result gives an expression for the condensation energy first 
derived by Wada [11]. 



m 

I - Jt% tl +$2 

+ \uo m y v s ^ <*2. ) (it) 

Let us clearly state what the Wada expression calculates. It is the strong coupling result 
for the condensation energy at thermodynamic equlibrium, under the assumption that there 
are no appreciable changes in the bosonic mode between the normal and superconducting 
states. In other words, it accounts for the appearance of the pairing vertex, as well as any 
changes to the fermionic self energy, but ignores any feedback effects between bosons and 
fermions. 

An equivalent expression for E c , more advantageous for numerical calculations due to a 
faster convergence at high frequencies was obtained by Bardeen and Stephen [10]. They 
noticed that there exists an integral relation between E^ and Es )U , 

N f nTj2 E^-pSd— - = (18) 

that in turn is the consequence of the fact that 

T W d 2 k^ m G s ^Jk) = T W d 2 k^ m G N ^ m {k) (19) 

m J m J 

as both quantities can be re-expressed as a cross-product 

T 2 £ / d 2 kD^ m G N ^ m {k)G s ^{k) (20) 

m,n 

Using equation (18), Bardeen and Stephen obtained the following expression for E ce i = E c . 

Ecei = -N f nTj2 (V% m + $L - |S SiW j) fl~ , f^' ) (21) 
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The practical importance of Bardeen-Stephen result is that at high frequencies, when 
^s,Lu m ~ ^N,ui m ~ Wm, the integrand in (21) behaves as $f, m /o;^, and the frequency sum- 
mation rapidly converges. This makes the Bardeen-Stephen expression more convenient for 
numerical computations than Wada's expression. 

We emphasize again that the equivalence of the two forms for E c is the consequence of 
the fact that in the phonon case, the change of the bosonic self-energy H(oo) between normal 
and superconducting states can be neglected, and the accuracy of this approximation is 
governed by the same parameter \v s /vp as the accuracy of the Eliashberg theory. In the 
next section we show that this is not the case for spin mediated pairing. We will see that for 
superconductors with an electronic pairing mechanism, the feedback on the pairing boson 
plays a crucial role, and the Wada's expression for E c would give competely erroneous results. 
Instead, the full expression 3 must be used. 

We pause now to connect with the BCS result for the condensation energy. BCS is a weak 
coupling theory. It assumes that the only change between the normal and superconducting 
states is the introduction of the pairing vertex $ which is given by the BCS gap A. (&bcs = 
A). The fermionic and bosonic self energies are both taken to be negligable. The BCS 
condensation energy is therefore calculated from the Wada equation with the substitution 
Suj m = u) m + S Wm PS uj m . Taking the zero temperature limit we obtain after some simple 
algebra 



This is the result that we already cited in the Introduction. We also see that the frequency 
integration in (22) is confined to u ~ A, i.e., the condensation energy comes from fermions 
in a narrow region around the Fermi surface. Although this result looks rather straightfor- 
ward, the issue of which fermions contribute to the condensation energy in the BCS case is 
nontrivial, and we discuss it in detail in Appendix B. 

III. MAGNETIC SUPERCONDUCTORS 

We now proceed to the case of magnetically mediated pairing. The bosonic mode that 
mediates the pairing is now the low-energy spin susceptibility. The Luttinger-Ward for- 
malism, which deals with an arbitrary bosonic mode is still valid, i.e., the Free energy has 




(22) 
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the same form as in Eq. 3. The only immediate modification is that now the second term 
of Eq. 3 has an extra factor of 3 reflecting the fact that all three components of the spin 
susceptibility contribute equally to the pairing. The Free energy then has the form 

n = -2Tj2{^og[el + t 2 (k) + ^ 2 (k)]-tJ:(k)G(k) + i^(k)F(k)} 
p 2 

+ 3 -Tj2{log[D-\ q )] + U( q )D(q)} 

Z q 

+ T 2 J2gLAG(k) X (k - k')G(k') + F(k) X (k - k')F(k')} + ... (23) 

k,k' 

The dimensionless bosonic propagator D(q) is now related to the magnetic susceptibility 
Xijio) = x{o)$ij as D(q) = x(<?)/x(Q) 0) where Q is the momentum at which the static 
susceptibility is peaked. Similarly to phonons, x(o) i s related to the bare susceptibility by 
X 1 ^) — Xo ~~ n(?)- We discuss the exact form of x(o) below. Finally, g q is now the 
coupling between the fermionic propagator and bosonic (magnetic) mode. 

At this point the formalism is quite general, taking into account all changes to both the 
fermions and bosons. In order to proceed further, we will need to assume a specific model 
that will allow us to neglect higher order terms in Eq. 23. We choose the spin-fermion 
model in which fermions are paired via their own collective spin excitations. Several authors 
have demonstrated that the exchange of collective spin fluctuations peaked at or near the 
antiferromagnetic momentum Q = (it, tt) yields an attraction in the d x 2_ y 2 pairing channel. 
We will be studing E c for this kind of pairing. 

The spin-fermion model is described by the effective action 



S = - I drf dr'J2cl(r)G^(k,T-T') C]ia (T') 

+\ jf dr £ dr> ^ Xo '(q. * ~ O S q (r) • S_ q (r') 

+g q [ P drJ2^(r).S^(r), (24) 

JO q 

where G^ik, r) = d T — v fe (k — k F ) is the bare Fermionic propagator, c^ a is the fermionic 
creation operator for an electron with crystal momentum k and spin projection a, s = 
c'ac, where <7j are the Pauli matrices, and g is the coupling constant which measures the 
strength of the interaction between fermionic spins and the collective spin degrees of freedom 
described by bosonic variables S q . For simplicity, below we assume g q — g to be momentum 
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independent. The bare spin susceptibility Xo 1 ^ 7 ") i s assumed to be peaked at Q and has 
a standard Ornstein-Zernike form, i.e., its Fourier transform over r is 

x(Q) 

Xo(q ' u) = i+aq-Q) 2 -(-/«-T ■ (25) 

where £ is the magnetic correlation length. This bare susceptibility comes from fermions 
with energies comparable to E F and should be considered as an input for the low-energy 
theory. 

The dimensionless coupling constant for the model of Eq. (24) (defined such that 
S(cj, kp) = \oj in D = 2) is, 

\ = 4u/(3v F C 1 ), (26) 

where u> = 9g 2 x(Q) / (167r£ 2 ). The numerical factors are choosen for further convenience. 
This overall scale oj and the coupling A are the only two parameters that matter at strong 
coupling. Other parameters, e.g., v s , turn out to be irrelevant (see below). Note that u is 
in fact independent of £ as x(Q) by itself scales as £ 2 . 

Near a magnetic transition, £ is large, i.e., A > 1, and spin-fermion model is a strong 
coupling theory in which feedback effects between fermions and bosons are extremely im- 
portant. It has been discussed in depth in [16] as a theoretical model as well as with respect 
to the cuprates, and provided an explaination of many unusual properties of the cuprates 
such as non- Fermi liquid behavior in the normal state [17], d x 2_ y 2 pairing [18], and the pseu- 
dogap [21]. We note, though, that there is still a great deal of contraversy regarding the full 
description of these phenomena. Since, however, it is generally believed that the cuprates 
are strongly coupled superconductors, it is instructive, regardless of one's prejudices, to 
examine the condensation energy for the spin-fermion model in detail, so as to illustrate 
the importance of properly accounting for all feedback effects when dealing with a strongly 
coupled system. We will show that the spin-fermion model accounts for many aspects of the 
experimentally measured condensation energy. 

A. The validity of the Eliashberg approximation 

We begin by briefly discussing the validity of the Eliashberg approximation for anti- 
ferromagnetically mediated superconductivity and how the assumptions inherent in the 
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Luttinger-Ward condensation energy formalism are justified. For the purposes of calcu- 
lating the condensation energy, there are two main issues to be discussed. First, as with 
phonons, it is possible in the spin-fermion model to separate the perturbative series such 
that the terms resulting from vertex corrections and higher order diagrams are " small" and 
therefore irrelevant. Second, since the effects of pairing are greatest near "hot spots" (points 
on the Fermi surface connected by Q = (tt,it)) one can invoke an effective momentum in- 
dependence for the problem, while still retaining the d x 2_ y 2 pairing symmetry. We briefly 
emumerate the reasoning leading to the above conclusions here. A more detailed discussion 
can be found in Ref.[16]. 

1. Spin fluctuations are collective modes of fermions, hence there is no difference between 
the Fermi velocity and the spin velocity, i.e., v s ~ vp. Then Xv s /vp ~ A, i.e., there is 
no way to separate perturbative series based on the difference between velocities. From 
this perspective, there is no Migdal theorem for spin fluctuations, and the perturbation 
theory with the bare spin propagator just holds in powers of the coupling A. 

2. The absence of small v s /vp implies in turn that the polarization operator U(Q,cu m ) = 
n Um is not negligible, as it is for phonons, but is rather dominant for A ^> 1. The 
consequence of this is that one must simultaneously solve for both the fermionic and 
bosonic self-energies. 

3. In the normal state, n Um oc u m at low frequencies, i.e., when IT Wm dominates the 
frequency dependence of the spin susceptibility (strong coupling), spin fluctuations 
become diffusive. This transmutation of the spin dynamics from propagating with 
v s ~ vp for A <C 1 to diffusive for A > 1 implies that at strong coupling bosons become 
soft compared to electrons. Such softness of bosons, is precisely the physics behind 
the Migdal theorem. Not surprisingly then, the diagramatic series for fermionic E(/c) 
obtained with a diffusive bosonic propagator again can be separated into two different 
subsets of terms: one set of terms now scales as powers of log A instead of powers of A, 
and the reduction of the expansion parameter is a direct consequence of the softness 
of bosons compared to fermions. As for phonons, the series in log A gives rise to vertex 
corrections and to the renormalization of the quasiparticle residue. There are also 
terms that form series in A. As with phonons, these terms come from boson-induced 
interactions between electrons and their own zero-sound modes. That these series 
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hold in powers of the same A as the perturbation series with a bare boson propagator 
can be easily understood as at low frequencies, the interaction between fermions and 
their zero-sound modes is mediated by a static boson, and hence is unsensitive to any 
transmutation of the bosonic dynamics. 

4. This separation of terms into perturbative series of A and log A allows an approximation 
similar to that made for phonons to be made here. In the magnetic case one neglects 
terms 0(log A) compared to terms of order A. This is not as good of a approximation as 
the neglect of \v s /vp terms for phonons as log A is also large when A is large. However, 
the numerical prefactors for log A series turn out to be small (a vertex correction is 
only (1/8) log A), and in practice the neglect of logarithmical terms is well justified for 
all physically reasonable A (A ~ 1 — 2 at optimal doping). In addition, from a purely 
theoretical standpoint the log A terms can be made parametrically small by intoducing 
a large number of fermionic flavors M (a vertex correction is then (l/8M)logA). 
Furthermore, a one-loop RG analysis of the logarithmical terms shows that they give 
rise to fractional exponents, but do not change the physics, and, in particular, do not 
affect the pairing problem. 

5. As with phonons, the series in A yields a T*(k,ou m ) that is predominantly dependant 
on frequency. More specifically, near A;— points on the Fermi surface connected by the 
antiferromagnetic wave vector Q (hot spots), T,(k,ij) depends on u, but not on e^. 
This momentum independence is crucial for the computation of the spin polarization 
operator: for E(fc,u;) = £(u;), the density of states is flat, and hLj m turns out to 
be independent of S(cu) and is the same in the normal state as it would be for free 
fermions: 

n( Wm ) = ^ = 4A 2 ^ (27) 

UJ s f UJ 

Here we have introduced the notation u s f = u)/(4A 2 ). This uo s f scales as £~ 2 and 
vanishes at the magnetic transition. 

6. Away from a hot spot, this independence from e& (i.e., on momentum perpendicular 
to the Fermi surface) prevails, but Yl(k,u m ) still depends on the momentum along the 
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Fermi surface. At T = 0, the self-energy takes the form 

E(k, uj) = X(k) 2UJ (28) 

where 

A(fc) = A/(l + (*0 2 ) 1/2 , cu sf (k)=u; sf (l + (k) 2 ). (29) 

and k is the component of k— along the Fermi surface. This k dependence cannotbe 
neglected at the lowest frequencies as near the transition as k appears in a combination 
with £. However, for uj ^> uj s f{k), the k-dependence disappears: H(k,u) ~ {iujuj) 1 / 2 
(we used the fact that 2\(k)(uj s f(k)) l l 2 = uj). Alternatively speaking, at uj > uj s f(k), 
the whole Fermi surface acts as one big hot spot. In this range, the Eliashberg theory 
becomes applicable for all momenta. 

7. We see that whether or not the k dependence of the self-energy can be neglected 
depends on what the relevant uj and k are. For the pairing problem, a detailed analysis 
shows that typical frequencies are of order uj, and typical k are of order ujjv F . Then 
typical uj s f(k) are of order uj, i.e., the momentum dependence along the Fermi surface 
introduces corrections 0(1). These corrections have been checked in [20] and found to 
be nonessential from the perspective of the basic physics. Note also that the theory 
assumes that uj < E F , otherwise the linearization of the dispersion near the Fermi 
surface would not work. This in turn implies that the pairing is confined to fermions 
in the near vicinity of a hot spot. 

8. In the phonon case, the momentum integration in the expressions for the Free energy, 
fermionic £ and anomalous vertex $ can be factorized and performed exactly. Such 
momentum related corrections are always small to the extent of Xv s /vf- In the mag- 
netic case, the corrections resulting from an analogous procedure are always smaller 
than 1, but whether or not they are small parametrically depends on the frequency. 
For frequencies relevant to the pairing, the corections to the factorization are again 
0(1). 

We see from the above considerations that at strong coupling, A > 1, the softness of 
fermions compared to bosons gives rise to an effective Migdal theorem, i.e., the vertex 
corrections are smaller than E which in turn predominantly depends on frequency. Contrary 

16 



to the phonon case, there is no single parameter governing the validity of the Eliashberg 
approximation. There are logarithmically divergent corrections, but they do not affect the 
physics of the pairing, at least in the one-loop approximation. There are also physically 
irrelevant 0(1) corrections stemming from the momentum dependence of the fermionic self- 
energy and the pairing vertex along the Fermi surface. An Eliashberg-type theory is valid 
when both corrections are neglected. As we stated previously, this is quite reasonable from 
a physical perspective, and we now proceed under the assumption that the momentum 
dependence of E and $ can be fully neglected. In the case of $, this implies that we 
approximate the d x 2_ y 2 pairing vertex (and, hence the gap A = + E(a;))) by its value 

at a hot spot, taking into account the fact that the d— wave symmetry of $ implies that it 
has a different sign between hot spots separated by Q. 



B. Thermodynamic potential at equilibrium 

We now proceed by calculating the thermodynamic potential at its equlibrium value. As 
in the phonon case, the condition that fl is stationary with respect to variations of E, $ and 
n gives 

= 3zTg 2 X (Q)J2G(k')x(k-k') 

k' 

= -3iTg 2 X (Q)^F(k')x(k-k f ) 
k' 

U(k) = -2Tg 2 X (Q) J2[G(k)G(k - k') + F(k)F(k - k')] (30) 

k 

Under the assumption of the momentum independence of E, $ and II the real and anomalous 
Greens functions have the form 

0„„(k) = fk " ^ 



e 2 k + EL + $ I 



1 



4 + + $L 

and the d x 2_ y 2 -pairing implies F w (k + Q) = — F w (k). Furthermore, as we discussed, the 
momentum integration in Eqs. 30 can be factorized and performed exactly. This yields 

E Wm = AvrT^ E - ' 
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n a/S 2 +$2 (1 " n u 



)l/2' 



4A 2 

n.JQ) = — TrT^I-i 

V V J_ <f> <f> 

_|_ ^^n + m 1 ^^n^^n + m j (32) 

a/S 2 + $ 2 a/1] 2 _|_ $2 ' 

We emphasize again that Eqs. (32) contain only two inputs: the overall energy scale a) that 
is set by the spin-fermion interaction, and the dimensionless spin-fermion coupling A oc £ 
that diverges as the system approaches the antiferromagnetic instability. We also recall that 
the energy scale u is the ultimate upper cutoff for the strong coupling behavior (S Wm < 0J m 
for uj m > u) m ), while dimensionless A can be represented as the ratio (2A) 2 = uj/uj s f of uj and 
another typical scale u s f that sets the upper boundary of the Fermi-liquid behavior in the 
normal state. We illustrate the form of S^ m , <J> Wm and II Wm in figure 1 for both the normal 
and superconducting state. 

Substituting equations 30 into 23 we obtain the equlibrium thermodynamic potential in 
a magnetically mediated superconductor: 

<Pk 
(2^ 

3^^ f d 2 k x {Q,o) n , n x{q,um) 

(2 7 r) 2|1Ogi X (g,c m ) J+iW X (Q,0) 

The difference between the above equation and equation 12 (apart from the extra factor of 
3) is that we have retained the bosonic part of the free energy as it is by no means small. 
In light of this, it is convenient to represent Q as the the sum of two parts Q e i comprising 
the "electronic" contributions and Vt spin comprising the "magnetic" part. Accordingly, 

E c = E Cte i + E CjSpin (34) 

where 



m 



^ = -TE / 7^{log^ + S 2 m + $ 2 J-^ m G Wm (A;)+z<|.^ m (A;) 



n spin = -T^ / ——{log — + Ii Wm } (35) 

The electronic term E CjC i = f2f z — Q% accounts explicitly for the appearance of the anomalous 
pairing vertex <J> Wn , and for the feedback changes to the fermionic self-energy. This term by 

18 




FIG. 1: Matsubara frequency solutions at A = 1 for S (a), <J? (b) and II (c) in the spin fermion 
model (eqs 32) for both the normal and superconducting states. Note that £ is a strong function 
of frequency and may not be neglected. The apparent non-convergence of £ at high frequencies 
is spurious and is discussed in the text. Further note that II changes appreciably between the 
normal and superconducting states. This change must be taken into account when calculating the 
condensation energy. 

itself leads to the Wada result for the condensation energy. The term E c spin = Q s S p in — Q™ pin 
accounts for changes to the spin propagator via the changes to the spin polarization operator 
Il Wm . Together these two expressions account the feedback effects between fermions and 
bosons in a strong coupling theory. We point out that the distinction between Q el and fl sp i n 
is quite artificial, as the two are intimately connected by mutual feedback. It is the sum 
of the two which is physically relevant, and the two parts of E c may not be considered 
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separately unless, as in the phonon case, one of them is negligable. 

As S and $ and II depend only on u, the momentum integration in Eqs.35 can be 
performed explicitly and yields 



+ \uj m \ l S ^"L s '" m + (36) 



E c>spin = U S,u; m - ^N,oj m + log " S,UJm (37) 

The first term E Cfi i is the Wada result. The second term E CySpin is new. Note that its 
expansion in II^ — II/v begins with the quadratic term. This is the obvious consequence of 
the fact that the Free energy is stationary with respect to a variation of II. 

C. A cancellation of divergencies 

At a first glance, the electronic part of the condensation energy is qualitatively the same 
as the phonon result. This turns out, however, not to be the case as E Cfi i in fact contains a 
divergent piece which is cancelled out by the divergence in E^ spin . Indeed, consider the high 
frequency part of E c ^i. At high frequencies, S Wm dominates over <& Wm) and the electronic 
part of the condensation energy reduces to 

Q el = -NfTrT^l^J - |E^J + ... (38) 

m 

where dots stand for other terms that are all finite, as one can easily demonstrate. Examine 
next the equation for £. By making the substitution A Wm = $ Wm o; m /S Wm we may write E Wm 
in the following form. 

= *w £ -7== (1 _ „ 1 w < 39 > 

Were the bosonic spectrum unchanged between the superconducting and normal states 
(Un = lis) then we could expand equation (39) in powers of A n , and would find that 
at large frequencies 

Z 8jUm ~ oc ^ (40) 

Since the gap A Wm is expected on physical grounds to vanish at the highest frequencies 
(and computations indeed confirm this), the frequency integral in (38) converges, i.e., the 
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electronic part of the condensation energy would be finite. The situation is very different 
when changes in II are taken into account. Although at high frequencies Hs,u n indeed 
converges to Hn,ui„, the two expressions are different at frequencies comparable to typical 
A Wm . Since for arbitrary large uj m in equation (39), there is a range of running uj n where II5 
and IItv differ, and E^r do not converge at high frequencies: ^s,u>„ remains larger than 
^n,uj„ by a constant. We illustrate this behavior in figure 1 (a). 

This non-convergence of E^v and E^ seems at first glance to imply an infinite result for 
the condensation energy! It turns out that this apparent infinity is compensated for by the 
spin part of the condensation energy. As written in equation (37) the spin condensation 
energy looks quite convergent if we use Tl^^ m oc u rn . However, the expression for the spin 
polarization operator is formally ultraviolet divergent, and extra care has to be taken in 
evaluating the difference between IIjv,uj m and Hs,w m - 

In what follows we explicitly re-express the divergent contribution in E Cfi i in terms of the 
spin polarization operator, and show that when we take the divergent piece from E ce i and 
add it to E c ^ S pi n , the dangerous Ii-N,um —Us,u> m term in E CjSpin is cancelled out, and the remain- 
ing terms are all convergent, and in evaluating them we can safely use the regularization in 
which the ultraviolent divergent piece in 11^^ is absent, and Ii-N,cj m oc u m . In practice, this 
reqularization amounts to evaluating the integral over first, and the frequency integral 
later. 

In order to accurately single out the divergent piece in E c ^ d and relate it to the spin 
polarization operator, we use a trick originally suggested by Bardeen and Stephen and define 
a mixed self energy E^s (NS stands for normal-superconducting). This is the normal state 
Eliashberg equation for E but with the superconducting polarization bubble. 



We plot Tjns given by the above equation along with Es in figure 2 and show that they 
converge at high frequencies. We then add and subtract E^vs from E Cfi i. 




n\T^2sgn(u n ) 



1 



(41) 



n 



(i - IWJV2 



in 
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FIG. 2: Spin-fermion solutions for S5 and T,^s as defined by equation 41 in Matsubara frequencies. 
Note that by replacing II at by II5 in the normal state expression for E, we have obtained a 
convergent expression for the electronic part of the condensation energy as discussed in the text. 

+ {|WJ-|^v,U} (42) 

One can easily make sure that E c>e i now consists of a convergent piece plus the divergent 
|Etvs| — |Sjv|- The above is actually \T, NS \ — |Ejv| as the a; m 's cancel. 

We now explicitly express this divergent piece in terms of the spin polarization operator. 
To accomplish this, we recall that the term |En| in E ce i arose from the integration over 
momentum of the EG term in Eq. 35. Writing this for both the normal and normal- 
superconducting self energies, we have: 

r d 2 k 

N f^ T Y^\^N,Ns,w m \ = T J2 tttt^^n^ns^Gn^ (43) 

m m J l Z7 v 

By using the expressions for E and II the above may be written as a term in E CySp i n as 
follows: 

d 2 k 2 r d 2 q 



T E I TTTTli-Zg^TY, f j^\xN,NS,u>MGN,u> n+m {k + q)}G N ^ n {k) 



3 f d 2 q 

= -7^— T J2 / TTT^X^NS^iq^N,^ (44) 
Performing the integration over q we find that 

-nTN f £ {\E NSjUm \ - lE^J) = ~T^- 2 £ 1^ log (45) 
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We now move the divergent piece from Q e i to Vt spin and write the condensation energy as 
E c = 5Q ei + 5Q spin with 



5Q el = -N f nT^(^h m + ^ m -\^NS,u 



+ \Wm\ ,~ n ) (46) 



3T ^„ „ „ _ i - n 



sn sptn = -— y: n 5 , Wm - n^ m + (i + n^j log ^ (47) 

The electronic part is now fully convergent. For the spin part, one can easily check that at 
large frequencies, when Hn,^ and Hs,uj m are both large, the expansion of the logarithm in 
E c ,spin cancels the dangerous Hs,u m — term. The remaining terms are all ultraviolet 

convergent, i.e., are insensitive to the regularization procedure used to evaluate IlN,u m - This 
implies that the condensation energy is actually free from divergencies, as it indeed should 
be based on physical reasoning. 



D. The relation between <5$7 e / and 5Q sp i n . 

At this point, the electronic and spin contributions to the condensation energy seem to 
be rather different as the electronic part contains Nf, while the spin part doesn't. However, 
5Q e i and SQ sp i n are in fact of the same order as we now demonstate. Indeed, as we already 
said, typical frequencies for the pairing are of order Cu, and at these frequencies 

n s ,^ ~ n N ^ ~ — (48) 

Similarly, 

^S,Un ~ ^N,U, n ~ oj (49) 

Then using equation 47 for 5Q, spin : 

S&spin ~ TK ~ -5- (50) 

where the last step uses the definition of u s f given previously. At the same time 

5Q el ~ NfU 2 (51) 
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The Fermionic density of states Nf is a product of 1/vf (the leftover of the integration over 
6^, and a typical k along the Fermi surface). As typical k ~ ou/vf, 



N f ~^ (52) 



Substituting this result into (51), we find 



uj 3 



5Q e i ~ — , (53) 
v F 

i.e., 5Q e i and 5Q sp i n are indeed of the same order. 

In the above discussion Nf appears as an extra parameter in 5Q e i. This is because in the 
calculations we neglected the momentum dependence along the Fermi surface (the actual 
momentum integral over dk is replaced by a typical k). If this momentum dependence was 
included (i.e., by using the self-energy from Eq. (28) with A;— dependent X(k) and u s f(k), 
then the electronic part would be free from uncertainties. Unfortunately, this computation 
also requires the knowledge of the k— dependence of $(&;) along the Fermi surface, which 
is technically difficult to obtain. In constrast, SQ sp i n is the result of a full two-dimensional 
integration over bosonic momenta, and the result for 5fl spin is free from uncertainties. 

Fortunately, it turns out that within the (approximate) computational scheme that we 
are using, Nf and l/(uo s f^ 2 ) can be related. Their relation follows from Eq. (45), as both 
the fermionic self-energy and the spin polarization operator are fully expressed in terms of 
lD and A. By evaluating the constant pieces in Sjv5 jWm — Sjv, Wm and in Hs,w m ~ ^-N,w m at 
high frequencies and comparing the two sides of Eq. (45), one can express Nf in terms of 
l/{ui s f^ 2 ). Once this is done, there is no further uncertainty in the condensation energy - it 
is given by the universal function of A times u) 2 /(^s/£ 2 ) oc u 3 /v F . 

At the risk of belaborating this point, we note that the fact that E c>e i and E c>spin are of 
the same order implies that one cannot replace Wada's expression for for electronic part of 
the condensation energy by a Bardeen-Steven formula. Indeed, the ability to do this in a 
phonon superconductors was based on the integral relation between £at jW and £s jW , Eq. 18. 
For our case, the r.h.s. of this relation is finite and is given by 

K = - -2/ du m \og J - —S{u m ) (54) 

where f^ m = Tl s ^ m /Tl n ^ m , and 

poo (f) 2 ruim V 

S(u m ) = / du n + / du n (55) 
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A simple order of magnitude estimate shows that K is of the same order as £ s , Wm , i.e., for 
magnetic superconductors, the transformation from Wada's to Bardeen and Steven's formula 
for the electronic condensation energy introduces a correction of the same order as E c ^. 

One final remark. The electronic and spin parts of E c are only of the same order of 
magnitude as long as Q < E F . When the effective coupling exceeds E F , typical k — 0(1), i.e., 
the whole Fermi surface is involved in the pairing. In this limit, the spin-fermion calculations 
are not controlable. Estimates show, however, that typical frequencies for the pairing now 
scale as uj s f£ 2 ~ v f /lu ~ J where J is the exchange integral for the corresponding Heisenberg 
model. (Recall that u> ~ g 2 x{Q)H 2 where the RPA approximation, g is equivalent to 
Hubbard U. In the same approximation, near a magnetic transitiion x(Q) ~ £, 2 /U, i.e., 
uj ~ U.) Estimating 5Vt spin at typical frequencies, we indeed find SQ spin ~ J in agreement 
with the result by Scalapino and White [4]. The same reasoning yields <55T2 e z ~ Nf J 2 ~ 
J 2 /v F <C J- We see therefore that at very large couplings, the spin part of the condensation 
energy clearly prevails over the electronic part, i.e., the condensation energy comes entirely 
from the spin part. This again agrees with Scalapino and White [4]. 

IV. THE COMPUTATIONS 

In this section we present our results for the electronic and spin contributions to the 
condensation energy for various A. In practice, we found it advantageous to perform the 
calculations of 5Q spin and Sfl e i in real frequencies rather than in Matsubara frequencies. The 
main reason for this was simply that we had previously evaluated £(u;), $(a;) and at 
real frequencies and various couplings and could use these results in the present computa- 
tions. A more subtle reason is that in retarded formalism, the problem of divergencies in 
E Cfi i and E c spin can be avoided in a straightforward manner, (see below) 

A. Condensation energy in real frequencies 

We first derive the expression for the condensation energy in real frequencies, in terms 
of retarded £(u;), $(u;) and IT(a;). The Matsubara equations for E Cfi[ and E CySpin given in 
equation 35 have the following form. 

E c , e i = -irTj2f( iuJ m ) 

m 
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E C;Spin = -TiT^g(iuj n ) (56) 

n 

where the E c ^ has a sum over fermionic frequencies and E c>spin has a sum over bosonic 
frequencies. The retarded form of these equations, assuming no branch cuts except on the 
real axis are 

E c , e i = - J /" ret (^)tanh(— )du 

f°° UJ 

E c , spin = - J f ret {u) coth(—)dw (57) 

where /" is the imaginary part of f (/' is the real part) and g is similar. It remains to 
analytically continue f(cu m ) and g(w n ) to the real axis. With the Matsubara definitions 
used in section 2, the analytic continuations are as follows. 

S(cj m ) -> -iS ret (u) 
$(w m ) -> <5> ret {uj) 

U{u m ) -> U ret (u) (58) 
The retarded formulas for the condensation energy are then 



E c , el = —Nf / {\p + ReE N (u)] 

J 



ImL s (uj) a - ReL s (u) f3 u 
+ -[1 ]}tanh-^ 

3 f 00 

+ /mlogi^M)}coth^ (59) 

Where \J <£> 2 (u;) — S|(w) = a + ifi. We point out that extreme care must be taken with these 
equations in order to get the correct sign of the imaginary parts of both the square root and 
the logarithm. 

We first point out that there is no divergent term in E cel . Indeed, in the Matsubara 
formalism, the divergent term comes from the fact that at high frequencies, Es(u; n ) and 
Ti N (u! n ) are separated by a constant. Since we defined S with an extra i, this constant is 
imaginary. On the other hand, the the first two terms in the retarded formula for E Cje i at 
high frequencies where $ — > can be written as 



ImyJ-Y? s (u) + ReE N (w) = ReZ N (uj) - i?eS 5 (cj) (60) 
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FIG. 3: The real part of the retarded self-energy for both the normal and superconducting states. 
Observe that S^ ret and H' N ret converge at high frequencies in contrast to the constant offset in 
Matsubara frequencies. The constant offset goes into the imaginary part of X re t which does not 
affect E c in the retarded formulism. 

This follows from the fact that £5(0;) = S'(cj) + z|£"(u;)| and the branch cut is on the 
negative real axis. We see that E c>e i only depends on the difference of i?e£ between the 
normal and superconducting states, and the integral of this difference is fully convergent. 
We illustrate this in figure 3 Analogous reasoning also shows that E cspin is also free from 
divergencies. 

Indeed, the absence of divergencies in the retarded formalism is just the consequence of 
using the Kramers-Kronig transform which misses the divergent pieces in E CtSpin and E Cj€ i. 
However, since we already demonstrated that the full 5Q is free from divergencies, we can 



safely use the Kramers-Kronig transformation separately for E c spin and E, 



The fact that no divergence exists for the retarded formulas also allows us to relate the 
prefactors in front of E c ^ e \ and E C)Sp i n in a straightforward manner. In real frequencies, Eq. 
45 takes the form 



f°° UJ 
—Nf I {_Re£jvs(u;) — ReT> N (uj)} tanh — duj 

J Zi J. 



7T 2 £ 2 JO 



{ReIlN(uj)Im log 



1 ~ n g (o;) 



1 — Hg(u) U) 

+ImlL N (uj)Re log . . } coth —duj 

1 — lljv(,Cc> J <£-L 



(61) 



Since in the retarded formalism, i?eE/vs(a;) — ReT lN (uj)duj is a convergent quantity, 
and the r.h.s. of (61) is also convergent, we can explicitly evaluate (numerically) both 
sides of (61) and relate Nf and 3/8n 2 uj s fC, 2 which we label as spin density of states iV s 
(N 8 = (8/3n 2 )u/v 2 F ). 
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FIG. 4: Electronic {E Cje i), spin (E CjSp i n ), and total (E c ) condensation energy per unit cell at T = 
for various couplings A. The lines are a guide for the eye. The sum of E c>e \ and E C;Sp i n produces 
a total condensation energy which is negative. We used Nf = Ist/eV, and N s ~ 0.17st/eV as 
explained in the text. 

B. Results 

As we have already stated, we use previously obtained results for £(a>), <&(a>) and n(a>). 
First, we computed both sides of Eq, (61) and evaluated the ratio Nf/N s for various A. We 
found that with very small variations, Nf/N s w 5.9. 

In Fig. 4 we present the results for the electronic and spin contributions to the conden- 
sation energy for different values of the coupling A. To set the overall scale, we adopt a 
commonly used estimate Nf = Ist/eV [19]. We emphasize that changing Nf will only 
change the overall scale and not the functional form of E C (X). 

We see from Fig. 4 that the total condensation energy is negative, as it indeed should 
be in a superconductor, but that this negativity is of a very different origin than in BCS 
theory. In BCS theory, which corresponds to A <C 1, i.e., u < u> s f, the system behaves as a 
conventional Fermi liquid. In this limit, the pairing potential is static, i.e., the spin part of 
E c is negligible, and condensation energy is entirely electronic and negative. We see however 
for A > 1, ie the strongly coupled regime, the the electronic contribution to the condensation 
energy is positive and quite large. From the figure it appears that the electronic contribution 
changes sign at A ~ 0.4. E c e i is negative below this coupling strength and is positive for all 
A > 1/2 presented in the figure. Second, for all A shown, the spin part E CtSpin is negative. 
It can be shown that E c ^ pin continues to be negative at A — > where II — > as A 2 , (see Eq. 
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32) Indeed, by expanding the logarithm in Eq. 61 we obtain 

uj c r°° 

E c , spin (\ -> 0) = -N s ^f / Im{U 2 - U 2 N } 

2 Jo 

= -N s iu sf / {ReUglmUs - ReU N ImIl N } < (62) 
Jo 

The above equation is negative, as in the retarded formalism, ReIl N = and both ReTl s < 
and Imlls < 0, and scale as A 4 as A — > 0. 

We also see from the figure that at large A, both the spin and the electronic parts of the 
condensation energy nearly saturate: to a large positive value for E c ^ x and a large negative 
value for E CjSpin . The total condensation energy E c is negative and much smaller than 
either E c>e i or E c>spin due to a substantial cancellation between these two components of E c . 
Although this cancelation seems quite delicate, it is actually robust since E Cfi i and E CjSpin 
are intimately linked via mutual feedback, and can not be considered separately, it is the 
sum of the two which has physical meaning. Any estimate of the total condensation energy 
based merely upon either the electronic or spin part will give a highly erroneous result. 

We now consider the functional dependence of E c on A. In figure 5, we see that the 
condensation energy flattens at A ~ 2, and decreases at large couplings despite the fact that 
the pairing gap increases monotonically with A [20]. This behavior is very counterintuitive 
from a BCS perspective, where the condensation energy scales with A 2 . It clearly indicates 
that for A > 1, the physics is qualitatively different from BCS theory. To emphasize this 
strong deviation from BCS theory we plot in Fig. 5 the strong coupling result of E c along 
with the BCS condensation energy —NfA 2 /2 using the same A and Nf. We clearly see that 
for A > 1, corresponding to optimally doped and underdoped cuprates, BCS theory yields 
qualitatively incorrect results for E c . 

Our results are in line with earlier work which demonstrated that for A > 1, the pairing 
predominantly involves fermions located in the non-Fermi liquid frequency range. For these 
fermions, retardation effects not included in BCS theory become dominant. Such retardation 
effects take place between the "upper" - uj and "lower" ou s f scales of spin-fermion theory. 
As u — 4X 2 u s f this ratio grows quickly, and already uj/u s f = 4 at A = 1. This explains 
why the deviations from BCS behavior are already strong at this coupling. Understanding 
in detail the strong coupling physics behind the decrease in E c is currently the subject 
of a separate study [21] and a complete theory of this phenomenon does not exist at the 
moment. Most likely, however, this decrease is a reflection of the fact that as A increases, the 
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FIG. 5: Total condensation energy E c compared with the BCS result E^ cs at T = for various 
couplings A. We used N s ~ 0.17si/eV and Nf/N s ~ 5.9 as explained in the text. Observe that 
the BCS condensation energy monotonically increases as the coupling gets larger, while the actual 
condensation energy flattens at A ~ 2 and slightly decreases at large couplings. 

actual attraction between fermions goes down, retardation of the spin-mediated interaction 
becomes the major factor, and the pairing process increasingly involves incoherent (diffusive) 
fermions and on-shell bosons. As the exchange of on-shell bosons is an energy conserving 
process, it can not lead to a gain in E c . Such behavior is very counterintuitive from a BCS 
perspective, where the pairing emerges due to an exchange of virtual, off-shell bosons, and 
the condensation energy scales with A 2 . 

One final comment. Although magnetically mediated superconductors are often com- 
pared to dirty superconductors, we point out that at T = 0, the physics of the two is 
already qualitatively different. Analogies between the two are often made due to the fact 
that thermal spin fluctuations scatter at finite momentum transfer but zero energy transfer 
and act in the same way as non- magnetic impurities [20]. However, in a dirty superconduc- 
tor with non-magnetic impurities, the condensation energy retains its BCS form despite of 
the fact that the superfluid stiffness is renormalized down [22]. Obviously, this is not what 
we found. 

V. KINETIC ENERGY 

As we stated in the introduction, several groups have argued [5, 6] that the condensation 
energy is driven by a gain in the kinetic energy which at strong coupling is negative (in 




30 



contrast to BCS theory) because of a strong "undressing" of fermions which bear a greater 
resemblance to free particles in the superconducting state than they do in the normal state. 

In this section we consider, within our model, the change in the kinetic energy when the 
system enters the superconducting state. The conventionally defined kinetic-energy for an 
interacting fermionic system is 

E km = 2Tj2 J -0^ e k G^Jk) (63) 

where G UJm (k) is the full fermionic Green's function that contains the self-energy. Integrating 
over momentum and subtracting the normal state result from E k i n in a superconductor we 
obtain 

5E kin = 2N f nTj2 \/%^ + i L - I I (64) 

m 

In the BCS limit, A < 1, <& Wn = A, £ = 0, E Wn = u n and 

6E™ S = 2N f nT ]T ^ + A^ - |u; m | (65) 

m 

which is obviously positive and furthermore depends logarithmically on the upper limit of the 
frequency integration, which is u s f in our case (we recall that in the BCS limit, u s f » Hi). 
At T = 0, we have 

5E™ S = N f A*]og^ (66) 

In the same BCS limit, the potential part of the condensation energy 5E^ S is also loga- 
rithmically divergent, and to a logarithmical accuracy cancels out SE^ S . The subleading 
terms do not cancel and yield E^ cs = —NfA 2 /2. 

We now consider finite A. As before, we perform the computations in real frequencies. 
The analytic continuation of equation 64 gives 

r°° i ~ ~ to 

5E kin = +2N f ^ {Im V$ 2 (^) - £|M + ReH N {u)} tanh — dcu (67) 

The result of this calculation at finite A are presented in Fig 6. At low couplings the 
kinetic energy is positive, as one naively expects. At larger A, however, the kinetic energy 
passes through a maximum at A ~ 2 and then becomes negative at large A. 

As we already mentioned in the Introduction, the sign of E c depends on the interplay 
between two competing effects: the effect of particle-hole mixing that increases E kin , and 
the change in the self-energy due to the "undressing" of fermions that lowers E kin . At weak 
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FIG. 6: Kinetic energy 5Eki n compared with total condensation energy E c at T = for various 
couplings A. The parameters are the same as in Fig. 2. The kinetic energy change is positive at 
low couplings, but negative at high coupling. 

coupling, the particle-hole mixing obviously dominates. The sign change between small and 
large A implies that at strong coupling the situation is reversed, and the lowering of E^ n 
via the change in the self-energy due to the "undressing" of fermions overcomes the effect of 
particle-hole mixing This behaviour is very similar to that obtained by Norman et. al. [5]. 
As the first term in E c>e i is equal to —5E kin /2 (see Eqs. (37 and 64), one can indeed argue 
that the condensation energy at large couplings is at least partly driven by the lowering of the 
kinetic energy. However, a comparison of Figs. 4 and 6 shows that this is just another way 
to interpret strong coupling effects that affect both the fermionic and bosonic propagators 
via mutual feedback. 

A simple explanation of why this is so is the following. In the superconducting state, the 
spin decay into fermions is forbidden at energies smaller than 2A. This simultaneously gives 
rise to two effects. First, the spin propagator develops the excitonic (resonance) peak at 
u res < 2A. The energy released by the creation of an exciton results in a gain in the magnetic 
part of the condensation energy. Secondly, the fermions cannot decay until their frequency 
exceeds A + u res (this is the magnetic analog of the Holstein effect). The elimination of 
fermionic scattering at low frequencies implies that the fermionic self-energy E(a>) in the 
superconducting state is reduced compared to that in the normal state. This lowers the 
kinetic energy. Obvously, the two effects (the gain in the magnetic part and the lowering of 
the kinetic energy) come from the same physics. 
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VI. CONCLUSIONS 



Our goal in this paper was to emphasize the importance of taking all contributions to 
the condensation energy into account when considering a strongly coupled superconductor. 
Specifically we considered the case of d x 2_ y 2 pairing mediated by the exchange of near- 
critical overdamped antiferromagnetic spin fluctuations. We demonstrated that although 
Eliashberg theory is valid for a strongly coupled magnetic superconductor, the reason for its 
validity is qualitatively different from that for phonon superconductors as the spin velocity 
and Fermi velocity are of the same order. Due to this fact, approximations appropriate for 
phonon superconductors are generally not valid in magnetic superconductors. Specifically, 
we demonstrated that the assumption that the bosonic polarization bubble can be neglected, 
which was rigorously justified by Bardeen and Stephen for phonon superconductors, breaks 
down for magnetically-mediated superconductors and makes the Wada and Bardeen-Stephen 
formalisms invalid. 

We obtained the full expression for the condensation energy within the spin-fermion 
model and showed that the spin and electronic parts of the condensation energy are of 
the same order Q 3 /vp and both depend only on the dimensionless coupling A. The BCS 
behavior is restored at A < 1. Even at moderate couplings, the condensation energy is 
highly non-BCS. The electronic contribution to the condensation energy is positive, while 
the spin part is negative and larger in magnitude than the electronic part which makes 
the full E c negative. As in the BCS limit the electronic condensation energy is negative 
and equal to — iV}A 2 /2, this implies that the electronic condensation energy changes sign 
at a rather small A. We found that at large A, both the spin and the electronic parts of 
the condensation energy nearly saturate. As a result, the full condensation energy flattens 
at A ~ 2, and decreases at large couplings despite the fact that the pairing gap increases 
monotonically with A [20]. This behavior is very counterintuitive from a BCS perspective, 
where the condensation energy scales with A 2 and is also inconsistent with the behavior of 
E c in dirty superconductors. This behavior results from the fact that there is a substantial 
cancellation between the spin and electronic parts of the condensation energy and the total 
E c is thus substantially smaller than either the spin or electronic parts. 

We argued that the reduction of E c at large coupling is likely the result of the fact that 
at strong coupling the pairing is predominantly due to an energy conserving exchange of 
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on-shell (real) bosons as opposed to the BCS theory in which the pairing is caused by the 
energy non-conserving exchange of virtual bosons. The reduction of E c at strong coupling 
also indicates that in this limit coherent superconductivity becomes fragile, and a large A 
only indicates that the system needs a finite energy to destroy spin singlets. 

Finally, we computed the kinetic energy and found that at strong coupling it is negative 
which indicates that at high couplings the change in the self-energy due to the "undressing" 
of fermions, which lowers E^ n overcomes the effect of particle-hole mixing which tends to 
increase E kin . This behavior has no analog for phonon superconductors. We argued that a 
negative Ekm is fully consistent with a positive E c>e i, and that, in principle, it is correct to 
argue that the condensation energy at large couplings is at least partly driven by the lowering 
of the kinetic energy. However, our results show that the lowering of E c may be thought of 
as being due to either the the lowering of the kinetic energy, or the interplay between the 
lowering of AQ spin and the increase of E c e i. Both explainations are valid interpretations of 
the strong coupling effects which affect the fermionic and bosonic propagators via mutual 
feedback. 

By taking into account all contributions to the self energy, and taking Nf ~ Ist/eV as 
typical for near optimally doped cuprates, we obtained a small value for E c of ~ 15K at 
optimal doping, (which in our model corresponds to A ~ 1.5 — 2 [16], larger A describe 
underdoped cuprates). This is rather remarkable as all typical energies in the problem are 
much higher, i.e u ~ 2.5 — 3 x 10 3 K [16]. This small value of E c is partly due to small 
prefactors, but is also the result of substantial cancellation between the spin and electronic 
contributions to E c . Note also that Nf ~ Ist/eV that we were using is fully consistent with 
lD ~ 2.5 — 3 x 10 3 K. Indeed, using Fermi surface averaged vp ~ 0.6meV [20] we obtain 
N s = (8/37i 2 )u/v 2 F ~ 0.15 - O.lQst/eV. Using then N f /N s ~ 5.9 obtained in the paper, we 
find Nf ~ Ist/eV, i.e., precisely the same value as we used 

Our E c ~ 15K is in good agreement with experiment. Loram et al [3] extracted E c 
0.12/c^Tc ~ 10K from the jump of the specific heat at T c . The change of the functional 
form of E c around A = 2 is also consistent with the experimental fact that E c changes its 
behavior from BCS-like to non-BCS around optimal doping. A decrease of E c at strong 
couplings (i.e., for underdoped cuprates) is also consistent with what Loram et al found in 
the specific heat experiments in the underdoped regime [3]. We caution, however, that the 
relation between E c and the amount of the jump in the specific heat at T c , from which the 
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experimental E c was extracted may be more complex than in the BCS theory which was 
used to extract E c from the data. This analysis is clearly called for. 
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D. Pines, E. Yuzbashyan, and especially M. Norman for useful discussions. This research was 
supported by NSF DMR-9979749 (A. Ch) and by DR Project 200153, and the Department 
of Energy, under contract W-7405-ENG-36. (R.H.) 



VII. APPENDIX A - A RUNNING COUPLING CONSTANT APPROACH 

The condensation energy may also be computed using a general formula for the ground 
state energy of the interacting electron system [8, 13, 23]: 

£ - £ »=-'7o tt/w g - we: "' (68) 

where E is the ground state energy of free electrons, and the Green's function and the self- 
energy are evaluated for the running coupling constant Ai. This formula is valid both in the 
normal and superconducting state, The "effective" self-energy S* m is related to G UJm (k) in 
the same way as a conventional self-energy, i.e., as G~^(k) = ?(w m + S* m ) — e&. In the normal 
state, ^n,u> m = ^vmi where E is a conventional self-energy, while in the superconducting 
state, 

= ^ m + (69) 

Eq. (68) is particularly suitable for the strong coupling computations in the normal state. 
Here we can use 

r d 2 k 

J -^2 G "™( k ) = -inN f signw m (70) 

Subsituting this result into (68) we reduce E — Eq to a single frequency integral. Using the 
T = normal state result for the spin-fermion model S niWm = 2Acj m /(l + (1 + \co m \/ '^sf) 1 ^ 2 ), 
and and introducing the sharp upper cutoff for the low-energy theory at uj max ~ E F , we 
obtain 

e-e^-nM^I ' jjTTWTIr (71) 

This expression is convenient for the analysis of the variation of the ground state energy 
with A. 
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The condensation energy, i.e., the energy difference between normal and superconducting 
state, is given by 

Ec = ~i [ ^ / ^ (G s ,„jm„ m - G n ^(k)K.J (72) 
Using GT,* = - G 1 G), we can rewrite (72) as 

E c = N J £ *± (G s ,„Jk) - G n ^ m {k)){iu m - e k ) (73) 

Performing the momentum integration in the Green's functions, we obtain 



i: - - ^Fg** ' <74) 

This looks very similar to the Wada result ( Eq. 37) but note the minus sign in front of the 
u m . This extra minus sign results in a negative total E c for all couplings. 

The condensation energy, Eqn. (73) can also be formally divided into kinetic and potential 
energy terms, but this division is subjective for interacting systems, and we will not discuss 
this. 



A. BCS limit 

In the BCS limit, E c reduces to 



A d\i A 2 f°° duo. 



^ = -N f I ^ Ai, / — n (75) 



u 2 



/o X 1 ^Jo 

where A is the gap value for the running coupling Ai. Using the BCS relation between A 
and the coupling constant 

1 = X l N f / — 



one can rewrite (75) as 
From (76), 

Integrating over Ai we obtain 



E c 



°° dill 

-H AS - (77) 



A Al = A x e N f x N f x i (78) 



E c = -\ N f A 2 (79) 
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This is indeed the same result as we obtained using the Eliashberg formula. 

Note in passing that our previous assertion that the separation of E c into a kinetic and 
potential energy is subjective is true even in the BCS limit, as the interaction is the source 
of the pairing. Indeed, earlier we computed SEj, in in the Luttinger-Ward formalism and 
found that in the BCS limit, the kinetic energy scales as A 2 logw max /A (see Eq. (66)). In 
the running coupling constant formalism, the kinetic energy difference SEk in extracted from 
(73) in the BCS limit reduces to 



Using (76) and performing the computations with the logarithmical accuracy we find 



such that the sum of the two yields the correct total condensation energy, Eq. (79). This 
once again demonstrates that only total E c is a physically meaningful quantity. 



We numerically computed E c for A = 1 by calcuating the integrand of eqn 74 for A = 0, 
0.5 and 1, numerically fitting these three data points and then integrating the resultant func- 
tion over the coupling constant. This resulted in a condensation energy of E c = —15.8K 
compared to the —11.5K calcuated via the Eliashberg approach. This agreement is very 
impressive given the highly approximate nature of the calculation. Unfortunately, the run- 
ning coupling constant formalism is very computationally intensive, hence our motiviation 
for using the Eliashberg approach in this paper. 

VIII. APPENDIX B - A PECULIARITY OF THE BCS APPROXIMATION 

In this Appendix, we elaborate on our earlier discussions of the result for the condensation 
energy in the BCS limit. In BCS theory, the pairing problem may be described by an effective 
quadratic Hamiltonian 




(80) 



5E^ S = N A 2 /2, i.e., contrary to Eq. (66) the newly defined kinetic energy does not de- 
pend logarithmically on the upper limit of frequency integration. Similarly, SE^ S = —N A 2 



B. A = 1 



H - H + £k %,a a k,a + TT9aA a k, a a -k,(i + a k,pa-k, a ) 



(81) 



k,a 
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where g a p is the antisymmetric matrix [23]. The condensation energy can then be straight- 
forwardly obtained by simply averaging both the normal and anomalous terms in (81) Ex- 
pressing the average products of the pairs of operators in terms of frequency integrals of 
normal and anomalous Freen functions, and taking care to avoid double counting of the 
anomalous term, we obtain the zero temperature result 

E C = (82) 

where 

r deduj u 2 — e 2 



* J n (e 2 + u 2 )(e 2 + u 2 + A 2 ) ^ 
We know that the correct result is / = 1. However, a naive integration treating both oj m 

and e in (83) on an equal footing, results in a vanishing integral. This vanishing is surely 

artificial, as the 2D integral over deduj is logarithmically divergent and therefore the result 

does depend on the order of limits of the integration. This situation is a typical example of 

an anomaly. 

The correct way is to perform the integration over u first as in the Hamiltonian approach 
to the pairing, the interaction is independent of frequency, and hence there is no cutoff in 
the frequency integration. The integration over energy, on the other hand, is obtained by 
simplifying the original integral over momenta and is is only valid as long as the density of 
states is a constant which is true up to an upper cutoff A ~ ep. Integrating first over u m in 
infinite limits, and then integrating over we obtain 

/•A/A / , T 2 \ 

7 = 2/ o (^1-2^^=) (84) 

One can easily make sure that the integral converges at large x and therefore does not 

depend on the upper limit as long as A < A. Setting the upper limit to infinity we easily 

obtain I — 1 as it indeed should be. 

It is also instructive to reproduce the correct result by integrating over first. A formal 

integration in the infinite limit yields an incorrect I — —1, However, integrating in (83) first 

over €k between —A and A, and then over cu, we obtain 

2 roo / i + 2a; 2 A* A*\ 

I = / dx\ ~^^= tan" 1 -^^= - 2:rtan~ 1 — (85) 

where A* = A/ A. If we formally set A* = oo, the integral will be convergent and yield 
an incorrect result, I — —1. However, keeping A* large but finite we find after changing 
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variables to y = x/A* 

r = _! + V) 2 f *» ft-" 1 1 - *■>- 1 ' (V(A ' )2r ') (86) 

7T Jo v y y J 

Expanding under tan -1 in 1/(A*) 2 and evaluating the remaining integral we find 

I = -l + - r^- = -l + 2 = l (87) 
n Jo y + 1 

as it indeed should be. 

The physical implication of this result is that in the BCS theory, the condensation energy 
can be equally viewed as coming from the energy levels near the Fermi surface, as implied in 
Eq. (84) where the integral is confined to x = 0(1), i.e., to uo ~ A, or as coming from very 
deep levels below the Fermi surface, as implicated in Eq. (86) where the integral is confined 
to y = 0(1), or to uj ~ A ~ Ep. Such an uncertainty is typical for an anomaly which can 
be equivalently viewed as coming from either low energies or high energies [24]. 

This peculiarity, however, is only present in the BCS limit, where the gap remains finite 
even at the largest frequencies (this is what causes logarithmical divergencies in the inte- 
grals). In the Eliashberg theory, the gap vanishes at infinite frequency, the integrals are 
convergent, and the condensation energy can only be viewed as coming from the levels near 
the Fermi surface. 



A. A relation between E c in Refs. [5] and [12] 

The peculiar nature of the BCS limit also explains the apparent sign difference between 
the expressions for the condensation energy in Refs. [5] and [12]. In Ref. [12], Scalapino 
derived the condensation energy by averaging the interaction electron-phonon Hamiltonian: 

r rl d hrlt,i P 2 

where the last term is twice the expectation value of the ion kinetic energy. He then used 
Chester's result [25] for the relation between the isotopic dependence of the upper critical 
field and the change in the ion kinetic energy between normal and superconducting states, 
and found that for the isotope exponent a — 1/2, the change in the ion kinetic energy is 
precisely minus twice the change in the electronic propagator. As a result, the condensation 
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energy turns out to be minus the difference between the first terms in (88) in a supercon- 
ductor and the normal state: 

E c =< H sc >-< H n >= -i J ^^{oo + e k ) (G n (k,cu) - G s {k,u)) (89) 

[Note that we define E c with the opposite sign compared to Refs. [5] and [12] - our E c is 
negative, while their E c is positive]. 

For a nonzero self-energy E = £(u;), the momentum and frequency integral in (89) 
is ultraviolet convergent because the pairing gap A = A (a;) vanishes at w — > oo, and 
the ordering of the momentum and frequency integration does not matter. Performing 
integration over frequency first, one ontains, after simple manipulations, Wada's formula, 
Eqn. (17). In the BCS limit of vanishing £(u;) and a constant A, this yields E c = —N A 2 /2, 
as we discused earlier. 

The authors of Ref [5] , on the other hand, assumed that the condensation energy is not 
due to phonons, and that a is nearly zero, as in the near optimally doped cuprates. They 
argued that in this situation, the condensation energy should be given solely by the first 
term in (88), i.e., 

/d d kduj r d d k f° 

(90) 

where A(k,cu) = (l/n)ImG(k,u) is the quasiparticle spectral function. This expression has 
opposite sign compared to Eqn (89). Still, the authors of Ref. [5] argued that their expression 
also reproduces the BCS result E c = —N A 2 /2. 

The consideration in the previous subsection of this appendix shows that both expressions 
are correct in the BCS limit (despite the fact that the one is minus the other!). The 
way in which the condensation energy was obtained in ref [5] implies that the momentum 
and frequency integral must be physically motivated, i.e., if the momentum integral over 
d d k ~ N dtk is extended to infinite limits, one has to integrate first over momentum and then 
over frequency. On the contrary, the Wada expression for the condensation energy implies 
that the momentum integral comes first, and the frequency integration comes second. As 
we already know, in the BCS limit, interchanging the order of the integration changes the 
result by a factor —1. This explains why the two apparently opposite results for E c actually 
yield the same condensation energy. 
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A comparison of the two results shows that that in the BCS limit, the change in the ion 
kinetic energy merely sets the proper regularization of the ultraviolet divergent integral for 
E c . If one regularizes the integral by restricting the momentum integration to |e&| < A ~ Ep 
(this is what we called the physically motivated regularization), the ion kinetic energy may 
be completely neglected. Reversing the ordering of the integration is equivalent to imposing 
a constraint on the frequency integral. This constraint can only come from the frequency 
dependence of the phonon propagator, i.e., from the kinetic energy of ions. Obviously in 
this situation, the ion kinetic energy cannot be neglected. 

Away from the BCS limit, the correct way to proceed in the general case is to use the 
Luttinger-Ward-Eliashberg expression for Free energy, Eq. (3). This expression is valid for 
arbitrary Yt(k,u) and it also includes the full feedback on bosons. Note that the expression 
for Q is more general than the Eliashberg theory for phonon superconductivity as Eqn (3) 
actually doesn't assume that £ depends only on frequency, and that one can factorize the 
momentum integration in the expressions for the self-energy and the pairing vertex, Eqs (5). 
The only approximation in the Eliashberg formula is the neglect of vertex corrections (which 
account for higher-order terms labled by dots in (3)). Still, it is of interest to understand 
which of the two expressions for E c , Eqn 89 or 90, one should use for the situation where 
the feedback effect on the bosonic propagator can, for one reason or another, be neglected. 
We argue that the answer depends on whether the self-energy predominantly depends on k 
or on uj. 

Indeed, in the BCS limit, the physical motivation for the ordering of the momentum and 
frequency integrations (first over frequency and then over e^) was associated with the fact 
that, for a constant gap, the ultraviolet divergence of the momentum and frequency integral 
for E c was cut by that the integration over cannot extend to |e&| > A ~ Ep. Suppose now 
that the Eliashberg theory is valid, i.e., £ « S(cu), and A = A(u). As we already said, this 
A(a>) vanishes above some characteristic frequency uo which decreases when fermion-boson 
coupling increases. Once this frequency becomes smaller than A, the cutoff of the ultraviolet 
divergence is provided by the frequency integral, rather than momentum integral. In this 
situation, the physically motivated ordering of the integrations should be integration over 
€k first (in infinite limits), and the integration over frequency afterwards. This ordering is 
implicit in Eq. 89. Not surprisingly, this expression yields Wada formula which in turn 
follows from the Eliashberg Free energy provided that £ = £(u;). If instead we used Eqn. 
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(90), we would obtain a rapid variation of E c once uj becomes smaller than A, and eventually 
Eqn (90) would yield the result opposite in sign to Eq (53). 

On the other hand, when the self-energy predominantly depends on k, the cut of the 
ultraviolet divergence is still provided by the momentum integral. In this situation, the ionic 
kinetic energy can be neglected, and Eqn. 90 should be used. This can be explicitly verified 
by comparing the Luttinger-Ward-Eliashberg formula with the bosonic term dropped, Eqn. 
(16), with Eqn. (90). Setting E(k,uj) = S(e fc ) and &(k,u) = $(e fc ) in (16) and integrating 
over frequency we obtain 



where S e = e+E(e). This expression is the analog of the Wada formula for the condensation 
energy (17). We now do the same in Eqn. (90). The frequency integration is straightforward 
and by proper evaluation of the arguments of the logarithms we obtain exactly the same 
expression as Eqn. (91). This proves our point that in theories where the feedback from 
the pairing on bosonic propagator can be neglected, Eq. (16) is valid for any Tl(k,u), while 
Eqs. (89) or (90) are valid when H(k,u) ~ S(c<j) and Yl(k,u) ~ S(/c), respectively. 
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